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Abstract. Optical spectral variability of quasars and BL Lac Obiects is com pared 
by means of the spectral variability parameter (5 (Trevese & Vagnetti 2002). Both 
kinds of objects change their spectral slopes a, becoming bluer when brighter, but 
BL Lac Objects have smaller /3 values and are clearly separated from quasars in 
the a — P plane. Models accounting for the origin of the variability are discussed 
for both classes of objects. 
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1. Introduction 

Variability of the spectral energy distri- 
bution (SED) of Active Galactic Nuclei 
(AGNs) is a powerful tool to investigate the 
role of the main emission processes in dif- 
ferent AGN classes, and the origin of their 
variations. The most common behavior in 
the optical band is that AGNs become 
bluer, i.e. their spectrum becomes harder, 
when brighter. This has been shown for 



ent ly shared by BL Lac O bjects, as shown 



and Seyferts (Cutri et al. 


1985 




Edelson, 


KroUk & Pike 1990|; Kinney et al. 


1991; 


Paltani & Courvoisier 1994) and for one 



complete sam ple, i.e. for the 4 2 PG QSOs 
monitored by Giveon et al. (1999) in B 



and R for 7 years. Evidence based on two 
epochs has been found als o for the faint 
QSO sample in the SA 57 ( Trevese, Kron, 
BunonelpOOll) . The same trend is appar 



by D'Amicis et al. 



( ^002 ), who present 5- 
year long B, V, R, I light curves for eight 
objects. A quantitative estimate of color 
variability is necessary to compare the ob- 
servational data with emission models and 
to compare different AGN classes. In a. pre- 
vious paper ( [Trevese k Vagnetti | p002|) , we 



introduced the spectral variability param- 
eter (3 = Aa/Alog/^, /jy being the spe- 
cific flux and a = 91og/i//c?logi/ the spec- 
tral slope. From the average a and /3 val- 
ues, it was possible to derive constraints 
on the variability mechanisms. The spectral 
variability parameter was then es t imat ed 
dVagnetti, Trevese fc Nesci | ^002]) for th e 



BL Lac Objects of D'Amicis et al. (2002) 



Observations and models can be compared 
in the a — (3 plane, which is reported in Fig. 
1 for both the PG QSOs and the BL Lac 
Objects. 
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2. Quasars 

Various models have been compared with 
the spectral variability of PG QSOs 



(Trevese & Vagnetti 2002) and are shown 
in Fig. 1. 

(i) The dot-dashed line represents the 
spectral variability due to small tempera- 
ture changes for a sequence of black bod- 
ies with different temperatures (increasing 
from left to right). 

(ii) We evaluate the effect of the 
host galaxy through numerical simulations 
based on templates of the QSO and host 
galaxy SEDs, derived from the atlas of nor 



mal QSO continuum spectra (Elvis et al 



1994). We added to the fixed host galaxy 



template SED the average QSO spectrum 
with a relative weight measured by the pa- 
rameter T] = log(L^/Lfj), where and 
I/f^ are the total H band luminosities of 
the QSO and the host galaxy respectively. 
Variability is represented by small changes 
A77, around each 77 value, with an ampli- 
tude corresponding to a r.m.s. variability 
as = 0.16 mag in the blue band. The re- 
sult is shown in Fig. 1 for —3 < 77 < 3 and 
for z = 0, represented by a thin, continuous 
line, and clearly shows that the effect of the 
host galaxy is not sufficient to account for 
the observed changes of the spectral slope. 

(iii) We considered the ac c retion disk 
model of ^icmiginowska et al. ( 1995| ), cor- 
responding to a Kerr metric and modi- 
fied black body SED, which depend on the 
black hole mass M, the accretion rate M 
and the inclination 9. A change of m = 
Mc^/Le {Le being the Eddington lumi- 
nosity) produces a variation of both lumi- 
nosity and the SED shape. The result is 
represented in Fig. 1 by a thick, dashed 
curve, for rh varying between 0.1 and 0.3, 
and for 9 = (face on disk). The spectral 
variations are clearly smaller, on average, 
than the observed ones. This means that a 
transition e.g. from a lower to a higher M 
regime implies a larger luminosity change 
for a given slope variation, respect to what 
is observed. 



(iv) Transient phenomena, like hot 
spots produced on the accretion disk by 



instability phenomena (Kawaguchi et al 



1998 ), instead of a transition to a new equi- 



librium state, may better explain the rel- 
atively large changes of the local spectral 
slope. We use a simple model based on the 
addition of a black body flare to the disk 
S ED, represent e d by the average QSO SED 



of Elvis et al. 



(|99j). The result is shown 

in Fig. 1 by the large filled squares, cor- 
responding to hot spots with T « 10^ K 
(upper), and T « 6 • 10^ K (lower). 

3. BL Lac Objects 

The continuum spectral energy distribution 
of blazars from radio frequencies to X and 
7-rays can be explained by a synchrotron 
emis sion plus inverse Compton s catter - 



ing (Sikora, Begelman, & Rees 1994) 



Variability can be produced by an inter- 
mittent channeling into the jet of t he en- 
ergy pr o duced by the central engine. 3pada 



si al. I (2001) have considered a detailed 
model where crossing of different shells 
of material, ejected with different veloci- 
ties, produce shocks which heat the elec- 
trons responsible for the synchrotron emis- 
sion. The resulting spectra are compared 
with multi-band, multi-epoch observations 
of 3C 279 from radio to 7 frequencies, show- 
ing a good agreement. In the case of the 
eight objects of our sample, B, V, R, I 
bands are sampling variability of the syn- 
chrotron component. This component can 
be roughly described by a broken power 
law characterized by the break (or peak 
in lyL^) frequency and the asymptotic 
spectral slopes ai end 02 at low and high 
f requency respectively, ai > —1, a? < — 1 
( Tavecchio, Maraschi fc Ghisellini 1998 ) 



We adopt (Vagnetti, Trevese fc Nesci 2002) 



the equivalent representation 

2ZyT) 



(1) 

for a stationary component of the SED and 
we add to it a second (variable) component 
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Fig. 1. The spectral variability parameter (3 versus the average spectral slope a for two 
samples of QSOs (dots) and BL Lac Objects (open circles). QSOs & BL Lac Objects 
appear clearly separated. Model predictions for both classes of objects are discussed in 
the text: black body (dot-dashed line), effect of the host galaxy (thin line), change of 
accretion rate (thick dashed line), hot spots (large filled squares), synchrotron emission 
for BL Lac Objects (vertical lines). 



with the same analytical form but different 
peak frequency Vpi and amplitude Lpi to 
produce spectral changes: 



L„ ^ Lp- H{u] ai, a2, Vp) + 

+Lp, ■ H{v;ai,a2,Vp') (2) 



The addition of the second component 
mimics the behavior of the synchrotron 
emission of model spectra when shell cross- 
ing occurs, producing an increment of emis- 
sion with Vpi > Up, due to newly accel- 
erated electrons. With such a representa- 
tion we can compute a and /? as a func- 
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tion of Vpi/vp^ for different values of Up 
and given values of ai, a2 and LpijLp. 
We ado pt typical values of the asymptotic 
slope s ( Tavecchio, Maraschi fc GhisellinT] 
1998 ) a\ = —0.5, a-i = —1.75, we as- 
signed to Vp different values in the range 
10^'' — 10^*^, we made Vpi/vp vary in the 
range « 1 — 10 and we adopted LpijLp 
corresponding to a magnitude change of 1 
mag r.m.s. The results are shown in Fig. 1, 
where for a we use the slope of the station- 
ary component. The three vertical lines are 
computed for logz/p — 14,14.5,14.8 from 
left to right respectively. The computed 
lines fall naturally in the region occupied 
by the data, namely it is possible to ac- 
count for the position of the objects in the 
a — (i plane with typical values of ai, a2, 
Vpi jvp and z/p, corresponding to th e overal 
SEP of the objects considered (see 



et al. 1998) 



Fossati 



4. Conclusions 

We show that the spectral variability pa- 
rameter /? is a powerful tool to discriminate 
between different models of the variability 
of AGNs. Hot spots on the disk, likely pro- 
duced by local instabilities, are able to ac- 
count for the observed spectral variability 
of QSOs. 

We show that BL Lacs clearly differ 
from QSOs in their a,/3 distribution. A 
simple model representing the variability of 
a synchrotron component can account for 
the observed a and /? values. 

In the framework of wide field variabil- 
ity studies, we stress that observations in 
at least two photometric bands, repeated 
on the same field at many epochs, would 
allow a detailed test of variability models, 
extending our knowledge of the emission 
processes in AGNs. 
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